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SUMMARY 

Results of an investigation to determine average outside surface 
heat-transfer coefficients for a cascade of symmetrical impulse tur- 
bine "blades are given. It has been theoretically shown that the 
effectiveness of both direct and indirect turbine-blade cooling 
methods depends to a great extent on the value of the heat-transfer 
coefficient between the gases and the blades. Because of the lack 
of data necessary to establish fundamental heat-transfer laws appli- 
cable to turbine blades, the NACA has started a comprehensive pro- 
gram that will meet these needs. As a part of this program, the 
heat transfer in a cascade of rim-cooled symmetrical impulse turbine 
blades was investigated. The blades were heated at the roots and 
cold air was exhausted past them. Adiabatic tests were first made 
to determine the thermal recovery factor for the blade, so that the 
heat-transfer coefficients could be based on the difference between 
the blade and effective gas temperatures. 

The recovery factor was found to vary only slightly (from 0.78 
to 0.89) with an exit Mach number ranging from 0.3 to 1.0 and was 
independent of the Reynolds number. The Nusselt number based on the 
effective gas temperature was independent of the Mach number In a 
range from 0.3 to 1.0. The results of the heat-transfer tests can 
be represented within ±10 percent by 

Nu = 0.14 (Re) i °* 68 (Pr) 1 / 3 

where Nu is Nusselt number. Re is Reynolds number, and Pr is 
Prandtl number. The characteristic dimension is arbitrarily taken 
as the blade perimeter divided by n and the inlet Reynolds 
number is used. This equation expresses equally well the results 
of investigations made by the General Electric Company on a cascade 
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of impulse blades five times the size of the NACA blades. No effect 
due to the difference in curvature of the two blades was evident. 

The results of the low-temperature runs are shown to be applicable 
to the prediction of the effectiveness of rim cooling at elevated 
temperatures for the blades used. 


INTRODUCTION 

Cooling methods for gas turbines can ordinarily be classified 
as either indirect or direct. With the indirect method a coolant 
is applied to the blade root or the tip, the rest of the blade 
being cooled by conduction; with the direct method, a coolant is 
forced through passages in the blade. Results of theoretical 
analyses made at the NACA Cleveland laboratory (references 1 to 3) 
to investigate both types of cooling indicate that the amount as 
well as the effectiveness of the cooling depends to a great extent 
on the convection heat-transfer coefficient between the hot gases 
and the blades. 

With indirect methods, the cooling effectiveness, as measured 
by the allowable increase in effective gas temperature, increases 
with a decrease in the beat-transfer coefficient. With direct 
cooling, the dependence of the cooling effectiveness on the heat- 
transfer coefficient is more complex, involving the ratio of the 
coefficient on the hot-gas side to that on the coolant side. How- 
ever, the same general trend as that for indirect cooling applies. 
It is then clear that accurate cooling calculations for either 
direct or indirect methods depend upon knowledge of the heat- 
transfer laws between the gases and the blades. 

Although a great deal of heat-transfer data is available on 
the flow of fluids past plates, cylinders, and airfoils, none has 
been published on turbine blades. The General Electric Company, 
however, ran tests on a turbine -blade cascade using room-air - 
and steam-heated hollow blades and obtained 

Nu = 0.14 (Re) 0,68 (Pr) 1 / 3 


where 

Nu Nusselt number, based on blade perimeter divided by Jt 
Re Reynolds number, based on blade perimeter divided by k 
P r Prandtl number 
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The temperature at which the air properties were taken was not given; 
however, the temperature range involved was so small that either the 
bulk or film temperature could be used without appreciably changing 
the results. It was recognized in this work that there might be an 
effect due to the reversal of heat flow, that is, heat transfer from 
the blade to the gas instead of from the gas to the blade, and a 
correction factor was applied. This correction, which was 
originally suggested for heating or cooling of fluids flowing through 
pipes, was estimated to have decreased Nu by only 3.4 percent. 

(The correction does not appear in the previous equation. ) Similar 
work has undoubtedly been performed in other countries, but reports 
of the results are unavailable at present* 

Because of the lack of correlated data necessary to establish 
fundamental heat-transfer laws applicable to turbine blades, a 
comprehensive program hss been initiated at the NACA Cleveland 
laboratory that will meet these needs (reference 4). The problem 
is being experimentally attacked in two ways, by studying the heat 
transfer (a) in static turbine-blade cascades, and (b) in full-scale 
gas turbines. Cascades allow elaborate instrumentation, easy 
control of the factors involved in the convection and radiation 
processes, and rapid investigations of many blade configurations. 

The extremely large buoyancy forces due to rotation, which tend to 
increase free convection, and possible flow disturbances caused by 
the blades passing stationary nozzles are absent in cascades. The 
magnitude of these effects, though they are probably small, must be 
determined from research with actual turbines. 

The over-all program of research on turbine cooling conducted 
at the Cleveland laboratory includes an investigation of the laws 
governing the heat transfer by convection between the gases 
surrounding a cascade of symmetrical impulse blades and the blades. 
The results obtained by using cold air as the fluid and by heating 
the blade roots are presented herein. Cold air, rather than hot 
gases, was used because its properties are well known, and there 
is no possibility of carbon deposits forming on the blades. By 
maintaining low blade temperatures, the influence of radiation was 
made negligible, and thus a fundamental relation for the convection 
heat transfer could be obtained. The blade temperature distribution 
was measured and the results were used in conjunction with a theory 
for the heat flow to compute the heat-transfer coefficient. 

The purpose of this investigation is (a) to present for the 
blades tested a relation between the Nusselt and Reynolds numbers 
describing the heat transfer over a variation of inlet Reynolds 
numbers (based on the blade perimeter divided by «) of from 
10,000 to 150,000 and a range of Mach numbers from 0.3 to 1.0, 

(bj to correlate and compare the results with those for other 
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turbine -blade cascades, cylinders, and streamlined bodies, and 

(c) to apply the results to predict the effectiveness of rim cooling 

by using the theory of reference 1. 


THEORY 

The general relation for forced-convection heat transfer is 
often written in the dimensionless form 

Nu = f (Re, Pr) = C (Re) r (Pr) 8 (l) 


where 

Nu 

h c 

d 

k g 

Re 

V g 

p g 

^g 

Pr 

c P,g 


V 

Nusselt number, -r — 
k g 

convection heat-transfer coefficient, (Btu/(hr) (sq. ft)(°F)) 

characteristic dimension of body, (ft) 

thermal conductivity of gas, (Btu/(hr) (ft) (°F) ) 

V CT dp c 

Reynolds number. 


v g UM g 
"6 


gas velocity, (ft/sec) 
gas density, (slugs/cu ft) 
gas viscosity, (slugs/ft-sec) 

Prandtl number, 

k g 

specific heat of gas at constant pressure, (Btu/(slug) (°F) ) 


The symbols used in the analysis are defined in appendix A. 


The Prandtl number can be considered a physical property for 
a given gas and it varies only slightly with temperature. Thus the 
geometrical shape does not affect it. For turbulent flow, the gas 

properties p g , c p,g> 811(1 k g that 83:0 needed for Pr > Re > 

and Nu are evaluated at the film temperature t f , which is 

assumed to be the mean of the average body temperature and the gas 

static temperature. The characteristic dimension d for a turbine 
blade is not immediately apparent and must be found by correlating 
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heat-transfer data for a number of different blades. The blade per- 
imeter or perimeter divided by it has been commonly used . The heat- 
transfer coefficient is defined as the rate of heat transfer per unit 
area per unit temperature difference between the gas and the blade 


The value of At depends upon the method by which it is evaluated, 
as will be indicated. 

The surface temperature is usually well defined and at low 
Mach numbers little difficulty is encountered in choosing a gas teim 
perature because the static and total temperatures are very nearly 
equal. However, at high Mach numbers, these temperatures differ 
appreciably and different coefficients can be obtained, depending 
upon which temperature is used. Regardless of which temperature 
(total or static) is used, h<; may be positive or negative for the 
same direction of heat flow and a value of zero must be assumed 
for h c for the adiabatic case. 

In order to overcome these disadvantages, the heat-transfer 
coefficient is based on an effective gas temperature t g ^ e , which 
is defined as the temperature a body assumes in the absence of 
heat transfer (the adiabatic body temperature). A coefficient 
based on this temperature is always positive and greater than zero. 
Furthermore, it has been shown to be independent of the Mach number 
and the temperature difference At. (See references 5 and 6.) 

For any particular geometrical configuration, the effective 
gas temperature can be related to the total and static temperatures 
by a recovery factor a, defined by the equation 


surface 



( 2 ) 


where 

Q rate of heat transfer, (Btu/hr) 

S surface area, (sq ft) 

At temperature difference between gas and blade surface, (°F) 


= ^.e - s = Ve - t g,.B 

' tg,t - *8,8 y//2Jc P(g 


( 3 ) 
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where 


t g>8 static temperature of gas, (°F) 
t„ + total temperature of gas, (°F) 

J mechanical equivalent of heat, 778 (ft-lh/Btu) 


In a theoretical study, Pohlhausen derived an expression for 
the recovery factor for flow past plates as a function of only the 
Prandtl number (reference 7). Numerous experiments have been per- 
formed with air flowing past plates, parallel to wires, normal to 
single cylinders, and inside cylindrical tubes to obtain recovery 
factors. (See bibliography listed an p. 4 of reference 5.) Results 
of these experiments indicate that there is only a small Reynolds 
number effect; the principal variable is the Mach number for a 
fixed Prandtl number. Eckert and Weise (reference 8) tested three 
shapes of turbine blades, for which they determined a. They found 
same variation with blade shape and Mach number, but made no mention 
of a Reynolds number effect. 

In equation (2), it is seen that once a is known for a par- 
ticular setup, At can be based an tg^ e by use 0 Q.ua tion (3), 
and it is only necessary to determine Q in order to find h c . 
However, it was difficult to measure Q accurately in the test 
apparatus used and an alternate method of computing h c was devised . 

The one-dimensional temperature distribution for a turbine blade 
heated or cooled at the root can be found from a heat balance of the 
blade (appendix B). If radiation is neglected, 


i cosh m (L - x) 
0 cosh mL 


(4) 


where 


6 


excess of blade temperature over effective gas temperature 
at blade position x, (°F) 


e, 


o 


excess of blade temperature at x = 0 (blade root) over 
effective gas temperature, (°F) 


m 
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L distance from blade root to point at which de/dx is 

zero, (ft) 

x radial position on blade measured from the root, (ft) 

h„ h convection heat-transfer coefficient from gas to blade, 

(Btu/(hr)(sq ft)(°F)) 

b blade perimeter, (ft) 

kjjj a average thermal conductivity of blade material, (Btu/(hr) (ft) (°F) ) 

A cross-sectional area of blade, (sq ft) 

The excess of the blade temperature over the effective gas tem- 
perature 9 is a function only of x, the quantities 0q, L, 
and m being constants for a given set of conditions. If the tem- 
perature distribution is known, that is, 9 as a function of x, 
then 6 q can be found by inspection, and m and L can be computed 
by use of the method of least squares. Because 



the heat-transfer coefficient can ultimately be determined. 

In deriving the expression for the temperature distribution 
(equation (4)), the following assumptions are made: 

1. The blade is of uniform cross-sectional area and perimeter 
over the blade height. 

2. The thermal conductivity of the blade metal, the average 
heat-transfer coefficient, and the effective gas temperature are 
constant over the blade height. 

3. The temperature gradients in any cross section of the blade 
perpendicular to the radius are negligible. 

4. Radiation can be neglected. 

5. At some position on the blade L, there is no heat flow 
and de/dx = 0. 

It will be shown later that these assumptions are actually 
experimental conditions in the setup used. 
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By using the experimental value of h c ^> Nu can be computed 
baaed on an arbitrary dimension for d. The exponent of Pr in 
equation (l) has been found to have a value of approximately one- 
third for the flow of a gas around a body. The constants C and 
can be determined by plotting Nu/(Pr)^/3 as a function of Re. 


DESCRIPTION OF APPARATUS 

The general arrangement of the apparatus Is shown in figure 1. 
Refrigerated air, under a positive pressure (5 in. Hg) , passed 
successively through a heater, a VDI orifice, a throttling valve, 
a calming section, an inlet nozzle, a test section, another throttling 
valve. And into the laboratory exhaust system. The air temperature 
could be held constant within ±0.5° F. 

Details of the test section are shown in figures 2 and 3. The 
blades are made of Inconel and are of a symmetrical impulse design 
with a constant cross-sectional area and perimeter to conform to 
the first assumption of the temperature-distribution derivation; 
they are brazed to a bronze dummy-wheel section with a pitch-line 
radius of 5.50 inches. The blade solidity, defined as the ratio of 
the blade chord to the pitch, is 1.92. The blades are shrouded at 
the tip to form a flow passage 1 inch in height. 

The blade assembly fits between a split Inconel nozzle block, 
only the exit half of which is shown in figure 3. The blade roots 
can be heated by conduction through the bronze wheel section. The 
heat is supplied by means of an electric furnace, the bronze pro- 
truding down into the furnace* 

The temperature distribution of one blade is measured by means 
of a radiation-type thermocouple probe shown in figure 4 and described 
in appendix B. A hole, l/l6 inch in diameter, extended radially from 
tip to root through the center of one blade. A l/l6-inch— inside- 
diameter tube, which acts as a guide for the probe, is fitted flush 
with the blade tip so that the two holes are concentric. By moving 
the probe in the blade, the temperature at any point can be determined. 
Probe positions are measured with a calibrated screw accurate 
within ±0.001 inch. Calibrations of the probe (appendix C) insured 
the accuracy of the readings. 

In order to check the validity of the assumption of a one- 
dimensional temperature distribution, fixed thermocouples were 
installed on three blades. Three thermocouples peened into 1/32-inch- 
diameter holes were in a plane perpendicular to the span of the blade 
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on each blade. The distances of the three planes from the tips are 
0.23, 0.63, and 1.20 inches. Chordwise locations of these thermo- 
couples are indicated by circles in figure 2. 

Air inlet and exit total temperatures were measured in 6-inch 
ducts immediately before and after the test section. These ducts 
are large enough, compared to the test section, that the total tem- 
perature can be read directly. The thermocouple probe and inlet 
total temperature were read differentially on a potent iometer in 
conjunction with a light-beam galvanometer. Inlet total pressure 
was measured by means of a probe located at the nozzle throat. 

Inlet static pressure was measured with wall taps. Exit total and 
static pressures were measured with two fixed calibrated tubes placed 
between two blades at their trailing edges, as shown in figures 1 
and 2. 


PROCEDURE 

Adiabatic runs to measure the blade recovery factor as well as 
heat-transfer runs were necessary in order to obtain heat -transfer 
coefficients based on the effective gas temperature. In both kinds 
of run, however, the same measurements were taken, namely , orifice 
conditions, blade inlet and exit pressures and temperatures, and 
the difference between the blade temperature and the inlet-gas 
total temperature as a function of the probe position. 

For the adiabatic tests, the air temperature was so adjusted 
that the blades assumed room temperature, and thus heat losses were 
minimized. The Reynolds number was held constant by holding the mass 
flow fixed, while the Mach number was varied from 0.3 to 1.0 by 
changing the pressure ratio across the test section. Runs were made 
with three different Reynolds numbers. 

For the heat-transfer runs, the blade roots were heated with 
an electric furnace as previously described, and unheated air was 
exhausted past the blades. Runs were made at: (a) fixed Reynolds 

number, in which the Mach number was varied as in the adiabatic runs, 
and (b) variable Reynolds and Mach number with a ten-fold range of 
Reynolds number. 


10 


NACA RM No. E8H12 


CALCULATION OF RESULTS 
Adiabatic Runs 


It is advantageous to write the recovery factor in the form 


a = 1 - ~ 


- t 


■'g^t 


g>8 


Y 2 
g 


Ea 


(5) 


2Jc 


P,g 


because the term (t^ ^ - t^ g ) was experimentally measured, and 

in the adiabatic run the effective gas temperature was the average 

blade temperature. The gas velocity V CT is a function of the total- 

& 

temperature and pressure ratio 


y g = 



r 2=11 


i 

1 — 

7-1 


( 6 ) 


where 

g acceleration due to gravity, 32.17 (ft/sec 2 ) 

y ratio of specific heat at constant pressure to that at 

constant volume for air, 1.395 

R gas constant for air, 53.30 (ft/°R) 

T_ a. total gas temperature, (°R) 

p static gas pressure, (in. Hg) 

P total gas pressure, (in. Hg) 

Because the total temperature was constant through the blade passage, 
either the inlet or exit velocity could be computed by equation (6), 
using the pressure ratio at that position. 

Results of total- and static-pressure surveys indicated that 
whereas the velocity distribution was uniform at the blade inlet, 
separation was occurring at some point in the passage j therefore 
a velocity gradient resulted from inlet to exit. The small size of 
the blades made it impossible to determine at which point the separa- 
tion took place. However, in analyzing the test data, it became 
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apparent that the conditions measured at the blade exit were con- 
trolling the thermal processes in the blades. Recovery factors were 
therefore computed using the exit velocity and correlated with the 
exit Mach number M Q 



(7) 


Heat-Transfer Runs 

The average heat-transfer coefficient has been related to the 
blade -temperature distribution by equation (4). In order to check 
the theory experimentally, the test setup must agree with the 
assumptions upon which equation (4) is developed. These conditions 
were met in the following manner: 

1. The blade was designed with a constant cross-section area 
and perimeter over its entire height. 

2. The temperature differences from blade tip to root were kept 
small enough that thermal-conductivity variations were negligible, 
less than 7 percent. A comparison of a more exact derivation in 
which the thermal conductivity was considered a function of the tem- 
perature and equation (4) where the average conductivity for the 
temperature range was used resulted in a negligible difference. Even 
at greater temperature differences than those used in the experiments, 
the error could be neglected . A constant gas temperature was achieved 
by providing ample mixing length in the large duct before the test 
section and by the use of unheated air . 

3. Fixed thermocouples in the blade, arranged to give a two- 
dimensional temperature distribution, indicated that gradients 
perpendicular to the radius could be neglected. 

4. In order to minimize radiation effects, maximum blada tem- 
peratures were kept below 200° F. Because of the characteristics of 
the temperature distribution, over half the blade was at a tempera- 
ture less than 100° F for all runs. The temperature of the blade 
surroundings varied between 50° and 75° F. 

5. A point of zero heat flow was obtained by maintaining the 
gas temperature below that of the room. In effect then, the blade 
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was heated at the root by the furnace and at the tip by the room 
air because it assumed approximately the gas temperature at this 
zero heat -flow position. 

The recovery factors found in the adiabatic runs were used to 
convert the measured temperature difference - t g ^ to one based 
on t g ^ e as follows: 

For a given M Q , found from equation (7), a value of a was 
chosen. The exit velocity was computed by use of equation (6) by 
using the total and static pressures measured at the blade exit. 

From equation (5), 

V ^ 

0 = (tm - ^e) = ^m " *g,t^ + (! “ a ) 2Jc p ^g 

where 

t m blade temperature, (°F) 

V„ - velocity of gas at exit, (ft /sec) 

With 0 known as a function of x, it was only necessary to 
determine the constants of equation (4) in order to determine h c . 
however, to apply the method of least squares to a nonlinear 
equation, the observation equation (4) had to be transformed into a 
linear equation by expanding in a Taylor's series, neglecting all 
terms of powers higher than one (reference 9). The resulting obser- 
vation equation took the form 



where 

m', L' corrections to approximate values of m and L found by 
least squares 

9 value of 0 computed from approximate m and L 

v 

For a given run, approximate values of m and L had to be 
assumed from which for each measured value of x was 

computed using equation (4). The results were substituted in equa- 
tion (8) and the normal equations were formed and solved for m' 
and L'. The corrections were added to the assumed values to give 
the m and L that best fit the data points. When either correction. 
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m' and L', was large, the process was repeated until a negligible 
change occurred. 

An approximate value of L was found by plotting the observed 
temperature distribution and picking the position at which d6/dx 
was equal to zero. The value of L varied from run to run. Because 


m 



b 

A 


and b and A are known, only h c ^t) and l% #a must be assumed 
to find m. The blade thermal conductivity, based on experimental 
data for Inconel, was chosen for the integrated average blade tem- 
perature. The average heat- transfer coefficient can be written 


*c,b 


Q 


f L 


b 

J i 


e dx 


-k, 


A|M 




r L 


d dx 




(9) 


where 


thermal conductivity of the blade material for the temper- 
ature at x=0, (Btu/ (hr) (ft) (°F) ) 



temperature gradient at the point x=0, (°F/ft) 


The numerator of equation (9) is found by measuring the slope 
of the observed temperature-distribution curve at the position x=0, 
and multiplying by -k^ 0 A. The denominator is obtained from a 
graphical integration multiplied by b. 


A typical temperature-distribution curve calculated from the 
theoretical equation obtained by least squares as compared to the 
observed data points is shown in figure 5. 


The total temperature at the blades was assumed equal to the 
average of temperatures measured immediately before and after the 
test section. Reynolds numbers based on both blade inlet and exit 
conditions were calculated. Inlet Reynolds numbers (Re) 1 were 
computed by dividing the orifice mass flow by the flow area. The 
characteristic blade dimension was arbitrarily taken as the blade 
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perimeter divided by n. Exit Reynolds numbers (Re) 0 were based 
on velocities and densities computed from conditions measured at the 
blade outlets. The exit static temperature was computed from the 
exit gas velocity using equation (6) and the exit pressure ratio. 
Gas thermal conductivity, viscosity, and Prandtl number were evalu- 
ated at the film temperature taken as the mean of the integrated 
average blade and static (inlet or outlet) temperatures. 


DISCUSSION OF RESULTS 

Recovery factor. - The recovery factor for the blade investi- 
gated is shown in figure 6 as a function of the exit Mach number. 

The average adiabatic blade temperature was obtained by integrating 
the local blade temperatures as measured by the thermocouple probe 
and dividing by the blade length. Three sets of variable Mach 
number runs were made, each at an approximately constant Re. The 
exit Mach number was varied from 0.3 to 1.0, and the average exit 
Reynolds numbers were 60,000, 120,000, and 132,000. The results of 
all the runs lie about a single curve and there does not appear to 
be an appreciable effect due to Re. For an exit Mach number range 
from 0.3 to 1.0, a changes only from 0.78 to 0.89. Below 
an M 0 of 0.3, the temperature difference (t g ^ " bg,e) became 
too small to measure with precision. In this range, the effective 
gas temperature can be assumed equal to the total gas temperature 
without seriously changing results. 

Also on figure 6 are shown some results taken from reference 8. 
Two impulse blades, designated A and B were tested in a cascade 
using air as the fluid. Blades A and B turned the air 128 
and 110°, respectively. Adiabatic blade temperatures were measured 
with thermocouples at the thickest portion of the blade section. 
However, it is not stated at what distance, say from the tip, the 
measurements were made. If the blade temperature were constant from 
tip to root, this omission would make no difference. However, 
additional tests showed that heat was being lost from the blade ends 
to the surrounding surfaces. The measured blade temperatures were 
therefore probably somewhat low and, in turn, reduced a for the 
blades. The Reynolds number range of the tests is not stated. Not- 
withstanding these possible discrepancies, reasonably good agreement 
was obtained between the tests reported herein and those of refer- 
ence 8. The greatest deviation, which is for blade A, is less 
than 10 percent. 

Pohlhausen (reference 7) has theoretically shown that a is 
a function of Pr and additional theory has indicated that this 
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function is /^Pr. For these reasons, a/ tf¥r is also plotted 
against M 0 in figure 6. Because all the data were obtained using 
air at room temperature, there is no relative displacement of the 
curves. 

Heat-transfer coefficients . - The results of the heat-transfer 
tests are shown in the conventional manner, that is Nu/(Pr)^/ 3 as a 
function of Re, in figure 7. The characteristic blade dimension 
is arbitrarily taken as the perimeter divided by it for Nu and Ee. 
The results can be written in equation form as 


with an accuracy of ±10 percent. An inlet Reynolds number range 
from 10,000 to 150,000 and an exit Mach number range Mq from 0.3 
to 1.0 are covered. 

The inlet Reynolds number is obtained by dividing the mass flow 
measured with the orifice by the flow area, and corresponds to con- 
ditions at the blade inlet. If the exit conditions are used to 
define an exit Reynolds number, the dashed curve results, which is 
based on the maximum gas velocity through the cascade and the 
corresponding density. The slopes of the two curves are approxi- 
mately equal, but Nu/(Pr)l/3 based on exit conditions is lower by 
13 percent and can be represented by 


with an accuracy within ±10 percent. 

It was analytically shown (reference 6) that the Nusselt number 
based on the effective gas temperature is independent of the Mach 
number. Verification of this theory is given in figure 8. Results 
of three constant (Re) 0 runs in which Mq was varied are shown. 

Inasmuch as it was impossible to hold (Re) 0 absolutely fixed, 
Nu/CPr) 1 / 3 was corrected for (Re) 0 variations from an average value 
for each run. The correction was applied by multiplying the 

observed Nu/tPr) 1 / 3 by the ratio of the average (Re) 0 to the 
observed (Re) 0 raised to the 0.63 power. In general, the experi- 
mental scatter for these runs is no greater than that shown in 
figure 7. In all three cases, the maximum deviation from the mean 
value is ±10 percent and only ±5 percent for most of the points. 

Effective gas temperature. - The use of heat-transfer coeffi- 
cients based on the effective gas temperature is a refinement needed 


Nu 


= 0.14 (Re)! 0 * 68 (Pr) 1 / 3 
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only when the temperature differences involved are small, as was 
the case for these tests. As a rather extreme example, consider a 
blade at 105° F, the gas total temperature at 80° F, and a Mq 
equal to 1.0. The ratio of h c ^ based on t g ^ ‘ to based 

on tg^ e is 1.64. If the blade temperature is increased to 1080° F, 
the ratio is reduced to 1.016. At higher gas -temperature levels, 
the general trend is the same, although the percentage differences 
increase . 


C0MPABIS0N OF EESULTS 

Turbine cascades. - The results of this investigation are 
compared with those of other forced-convection investigations in 
figure 9. The cascade data of the General Electric Company coin- 
cide perfectly with that of the MCA, based on (Ee)j^. Because 
General Electric based Ee on the mass flow measured by an orifice, 
the only proper basis of comparison would be (Ee)i, which is 
obtained in a similar manner. Both blades were impulse; the blade 
used by General Electric turned the air 134° and had a 3.08-inch 
chord; the MCA blade turned the air 106° and had a 0.68- inch chord. 
The solidity of the two cascades was the same. General Electric 
Company also assumed an average recovery factor a equal to 0.85 
for all runs. 

Notwithstanding the differences in curvature and size of the 
two blades, an excellent correlation was obtained by using the 
blade perimeter divided by « as the characteristic dimension. 

Of course, the blade perimeter or any constant times the perimeter 
could have been used equally well. No effect due to the amount of 
blade curvature was present insofar as these two investigations 
reveal. It is quite probable that separation occurred at about 
the same point in each cascade and thus similar flow was established . 
By investigating blades of smaller curvature than those used in 
these runs, that is, blades in which separation is delayed or 
eliminated, some effects of the curvature might be found. 

When differences in shape involve more than merely the amount 
of curvature, such as a reaction blade compared with an impulse 
blade, it seems unlikely that data will be correlated by only a 
blade dimension. 

Because convection heat transfer is a boundary- layer phenomenon, 
there is some likelihood that data of different blades will never be 
correlated by the simple means outlined. Possibly, only through 
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theoretical boundary- layer studies that take into account the veloc- 
ity distribution about the blade will methods be developed to calcu- 
late heat-transfer coefficients for an arbitrary shape. Another 
premising mode of attack is by use of Reynolds' analogy, which considers 
the similarity between heat transfer and momentum transfer. Whereas 
neither of these approaches is likely to be simple, they will be 
much more satisfactory than experimental determinations for every 
blade shape. 

Cylinders and streamline bodies. - The data for air flowing 
normal to single cylinders and streamline bodies (reference 10) are 
also shown in figure 9. Again, the slopes for all the curves are 
approximately equal, but the magnitudes of Nu/(Pr)-*-/ 3 differ. Values 
for single cylinders and those for streamline bodies are 20 percent 
and 10 percent lower than those for the NACA blades based on (Re)i> 
respectively. 

Application to rim-cooling theory. - In addition to the cold- 
air tests, a high-temperature run with hot gases was made. The 
results of the cold-air runs are used to predict the cooling 
effectiveness of the NACA blades for this run. 

The effectiveness of rim cooling defined in reference 1 is 
measured by the allowable increase in effective gas temperature as 
limited by the blade stresses. For a fixed blade life, speed, and 
amount of cooling, the effectiveness is a function of only the 
quantity mL. The value of mL will now be predicted from the flow 
conditions and compared with the test value for a high-temperature 
run. 


The experimental conditions are: 


Excess of blade temperature at root over effective 

gas temperature, 

Inlet Reynolds number • • • • • • ' 

Average blade thermal conductivity, Btu/(hr) (ft) (^F) . 

Blade length, ft • 

Effective gas temperature, °F 

Film temperature, °F . . . 

Gas thermal conductivity, Btu/(hr) (ft) (°F) 


Prandtl number 

Blade perimeter divided by Jt, ft 


-452 

49,600 

19.7 

0.1042 

1167 

1000 

0.036 

0.65 

0.0473 
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From figure 7, 


Nu 

(Pr) 


173 = 211 


h c ^ = 139 (Btu/ (hr) (sq ft)(°F)) 


and 


m = 43.4 ft" 1 

The geometrical blade length Is taken as equal to L. The pre- 
dicted value of mL is 4.53. From these values of m and L, 
the temperature distribution can be calculated using equation (4). 
The result is 


9 = -9.73 cosh 43.4 (0.1042 - 1 ) 

By using the experimental temperature distribution, the measured 
value of mL can be obtained by applying the method of least 
squares, as previously discussed. The value of mL thus obtained 
is 4.26. 

The allowable increase in effective gas temperature for the 
two values of mT. can be found from a cross plot of figure 6 in 
reference 1 . For a constant value of 9 q equal to -452® F , the 
following results are obtained: 



mL 

At g,e 



(<*) 

Exper iment al 

4.26 

190 

Predicted 

4.53 

183 


The predicted value of At g ^ e is 4 percent lower than the 
experimental value. Radiation effects are unaccounted for in the 
predicted case. Also L is taken as the geometrical length of the 
blade, whereas actually for this case it is slightly greater. In 
spite of these approximations, good agreement is obtained. In fig- 
ure 10, the predicted blade temperature distribution is compared 
with the experimental results. The maximum error is only 10 F at 
a blade temperature of 1047° F or 1.0 percent. The values of x 
and L in figures 5 and 10 are given in inches to simplify the 

curves . 
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SUMMARY OF RESULTS 

From an Investigation to determine average outside surface heat- 
transfer coefficients for a cascade of symmetrical impulse turbine 
blades, the following results were obtained: 

1. For a range of inlet Reynolds numbers from 10,000 to 150,000 
and exit Mach numbers of 0.3 to 1.0, the results of the heat-transfer 
tests can be represented within ±10 percent by 

Nu = 0.14 (Re)! 0 * 68 (Pr) 1 / 3 

where Nu is Nusselt number, (Re)* is Reynolds number at the 
inlet, Pr is Prandtl number. The blade perimeter divided 

by « was used as the characteristic dimension in Nu and Re and 
the gas properties were evaluated at the average film temperature. 

The heat-transfer coefficient was based on the difference between 
the blade temperature and the effective gas temperature. 

2. The preceding equation expresses equally well the result of 
cascade tests made by the General Electric Company on impulse blades 
five times the size of the NACA blade. No effect due to the differ- 
ence in curvature of the two blades was apparent. 

3. The effective gas temperature can be related to the total 
and static gas temperatures by a recovery factor, defined as the 
ratio of the difference between the effective and static gas tem- 
peratures to the difference between the total and static gas tem- 
peratures. It varied only slightly (from 0.78 to 0.89) with an exit 
Mach number ranging from 0.3 to 1.0 for the blade used and was 
independent of the Reynolds number. 

4. The Nusselt number based on the difference between the blade 
Anri the effective gas temperature is independent of the Mach number 
up to a value of 1.0, at least. 

CONCLUSIONS 

From the preceding results, the following conclusions can be 
drawn : 

1. The results of these low-temperature tests can be used to 
compute the effectiveness of rim cooling at elevated gas tempera- 
tures for the blades investigated. 
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2. Although a correlation was obtained "between two "blade cas- 
cades of different sizes "by use of the "blade perimeter divided "by 
« as the characteristic dimension in the Reynolds and Nusselt 
numbers, it is not clear that this dimension will hold for all 
blade shapes and configurations. Additional experiments are 
necessary to add to the limited amount of experimental data now 
available. 


Lewis Flight Propulsion Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland , Ohio. 
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APPENDIX A 
SYMBOLS 

The following symbols are used in this report: 
area of blade cross section, sq ft 
speed of sound, ft/sec 
perimeter of thermocouple plug, ft 
perimeter of blade section, ft 


Ci, Cg 

_ „ arbitrary constants 

c 3> c 4 


C, r, s 


c 


P 


F 

g 


h c 

h c,b 


J 

k 

L 

M 

m 


constants 

specific heat at constant pressure, Btu/(slug) (°F) 

characteristic dimension in Reynolds and Nusselt numbers, ft 

shape factor for radiation 

acceleration due to gravity, 32.17 ft/sec^ 

convection heat-transfer coefficient, Btu/(hr)(sq ft)(°F) 


convection heat-transfer coefficient from gas to blade, 
Btu/ (hr) (sq ft)(°F) 


mechanical equivalent of heat, 778 ft-lb/Bt 
thermal conductivity, Btu/ (hr) (ft) (°F) 
distance from blade root to point at which 
Mach number, V g /a 


rim- cooling parameter, 




h c,b b 
*k,a A ' 


ft 


-1 


u 


d0/d 


I 


is zero, ft 


Nu Nusselt number, h c d/k g 

p total pressure, in. Hg 
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static pressure, in. Hg 
Prandtl number, — ^ 



rate of heat transfer, Btu/hr 

gas constant for air, 53.30 ft-lb/(lbW°R) 


Reynolds number, 


surface area, sq. ft 
temperature, °R 
temperature , °F 


Telocity, ft/sec 

radial position on blade measured from root, ft 

thermal recovery factor relating effective, static, and 
total gas temperatures 


ratio of specific heat of air at constant pressure to that 
at constant volume, 1.395 


allowable increase in effective gas temperature due to rim 
cooling, °F 

excess blade temperature over effective gas temperature at 
blade position x, - t g ^ e ), °F 

viscosity, slugs/ft-sec 

density, slugs/cu ft 

Stefan-Boltzmann constant, Btu/(hr)(sq. ft)(°R) 



temperature difference between gas and blade surface, °F 


COSH DLL. 
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Subscripts : 
av average 

b blade 

c convection 

e effective 

f film 

g gas 

i inlet 

l value calculated for least-squares computations 

m blade metal 

o exit 

p ping of thermocouple probe 

r radiation 

8 static 

t total 

0 blade root (x=0) 

1 upper surface of probe plug 

2 lower surface of probe plug 

Primed symbols indicate corrections to assumed values found by 
least squares. 
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APPENDIX B 

DERIVATION OF ONE-DIMENSIONAL TEMPERATURE DISTRIBUTION 

A steady- state heat balance for a differential element dx 
between two planes perpendicular to the blade length and at a 
distance x and x + dx from the blade root will be assumed. The 
difference between the amount of heat entering the element by con- 
duction dQ^. and the amount leaving by conduction dQ( x + dx ) is 
equal to the heat leaving by convection dQ c if radiation can be 
neglected : 

dQ x - dQ( x + dx) = dQc 



If 


kc,b 

^m^a A 


m 


2 


and (t m - t g e ) is set equal to 0 , then because t g>e is 
constant, 

d 2 t m _ d£e 
dx 2 dx 2 


Equation (B2) becomes 


d 2 8 

dx 2 


m 2 0 


A solution of this equation is 


0 = cosh m (Cg - x) 


(B3) 


It is assumed that at soma position on 
(For most actual turbine installations, 
geometrical blade length.) At x = 0, 


the blade L, de/dx = 0. 

L is closely equal to the 

e = 9 q. 


Substitution of these boundary conditions in equation (B3) gives 
the final equation for the temperature distribution when radiation s 

neglected 

c osh m (L - x) ( 34 ) 


0 = 0 


0 


cosh raL 
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APPENDIX C 


THERMOCOUPLE PROBE 


The radiation- type thermocouple probe used to obtain the blade 
temperature distribution is shown in figure 4. A junction is formed 
by spot-welding 36-gage chromel-alumel wire to a small Inconel plug 
0.055 inch in diameter and 0.032 inch high. The wires are insulated 
by two-hole Alundum tubing that has a 0.032-inch outside diameter 
and extends to within l/2 inch of the Inconel plug. The gap is 
necessary to reduce the transfer of heat by conduction between the 
plug and the tubing. For mechanical strength, an Inconel tube is 
cemented around the upper portion of the Alundum tubing. 

Before experimentally calibrating the probe, the effects of 
plug dimensions and material, that is, plug height and thermal con- 
ductivity, on the indicated temperature were theoretically investi- 
gated. Because the probe fits snugly into only a l/l6-inch-diameter 
hole, it was assumed that all the heat transfer between the walls and 
the plug and wires occurred by radiation. 

The derivation of the temperature distribution in the thermo- 
couple probe plug will now be indicated. (See fig. 11.) A steady- 
state heat balance for a differential element dx on the plug 
between two planes perpendicular to the blade length at a distance x 
and x + dx from the blade root is assumed. The difference between 
the amount of heat entering the element by conduction dQ x and the 
amount leaving by conduction dQ^ x + is equal to the heat 

leaving by radiation dQ^.; if convection can be neglected , 


dQ x " dQ (x + dx) " dQ r 



= oBF (T m 4 - T p 4 ) dx 


(Cl) 


By expanding T 4 in a 
T a and by using only 


in a Taylor's series about an average temperature 
only the first two terms, the result is 


T 4 = 4T a 3 T - 3T a 4 


(C2) 


The wall temperature T m can be written 

_ m . cosh m(L-x) 

T m = T g,e + e 0 cosh mL 
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Equation (Cl) becomes 


d 2 T r 


*&> dx 2 


= aBF 


4T 3 T + 4T 3 cosh 
4T a T g,e a y 0 cosh mL 


4T a 3 T p 


(C3) 


d 2 T 


P = P 2 


dx^ 


T g,e + °P c ° 3h m ( L " x ) " T p 


(C4) 


A solution of equation (C4) is 


T P " T g> 6 


+ C 3 e 


px 


+ Ca© 


-Px 


p 2 <p 

m 2 -P 2 


cosh m(L-x) 


(C5) 


The constants C 3 and C 4 can be evaluated from the following 
boundary conditions: 

(a) At the probe position 


x = x x 

Measured at the top surface of the plug 

t p - t p,i 

where T„ n is the reading of the probe. 
P J X 

(b) Also at 

x = x-^ 


- -Vp (j£) 


x=xi 


The amount of heat radiated to the plug surface at xq is computed 
bv dividing the wall into layers and summing the total heat transfer, 
the radiation-shape factor for each position being taken n o 
arco.mt To this value is added the heat conducted through the 
thermocouple wires to give This amount of heat can he computed 

by differentiating equation (C5) using the proper cons an s o 
wire. 
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The reading of the probe Tp^i is assumed for a given position 
and blade conditions. As outlined, all the constants of equa- 
tion (C5) can be evaluated, and the heat transfer at the bottom of 
the plug at xg can then be calculated 


The first term Qg can be calculated in the same manner as 
omitting heat conduction through the wires. Values of ^p,l are 
assumed until equation (C6) is satisfied. 

From this type of calculation, it was concluded that: (a) The 

temperature drop through the plug was so small (less than 1° F) 
that the conductivity was of no importance and almost any material 
could be used; and (b) As the height of the plug was increased, the 
probe error increased. An optimum height of 0.032 inch was 
indicated. 

With the results of the analysis serving as a guide, several 
probes were built and calibrated, as follows: A l/l6-inch-diameter 

hole was drilled from tip to root through the center of a blade. 
Five holes, each l/32 inch in diameter, were drilled perpendicular 
to the original hole through the blade at different distances from 
the root. Thermocouples, 36-gage chrome 1-alumel, were inserted in 
the small holes so that they indicated the temperature at the 
surface of the l/l6-inch hole. A tube was soldered to the blade 
root through which cooling water was passed. The entire apparatus 
was placed in an electric furnace, and the readings of the probe 
compared with those of the fixed thermocouples. By varying the 
heat input, different gas temperatures could be simulated. 

It was found that probes in which the plug height was greater 
than 0.032 inch read lower than the true temperature, as was 
indicated by the analysis. By using a plug height of 0.032 inch, 
very satisfactory agreement was obtained between the probe and the 
fixed thermocouples. The probe was very sensitive and attained 
thermal equilibrium quickly. 



(C6) 
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Figure 1. - Experimental setup for rim-cooled turbine-blade cascade. 
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Hydraulic radius, in. 

0.092 

C ho rd , in. 

0.680 

Perimeter, in. 

1.785 

Cross-sectional area, sq in. 

0. 0797 

Pitch-line radius, in. 

5.50 

Humber of blades in 360° 

84 

Sol i d i ty 

1 .92 


Figure 2. - Cross sections of symmetrical impulse 

turbine blades. 
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Thermoc ouple - pr ob • 

guide 


Bronze dummy 
wheel section 


Cascade 
of b lade a 


NACA 


C- 14757 
4 - 22-46 


■pr* whj*- r ^=7 ; : , 

I Ilf l N 



Ft c / 


| Them mlo I * \ ^ 

■= Inconel exit- 1 

I seal ■•mVVijj 

block | 


(a) Unassembled. 



(b) Assembled. 

Figure 3. - Special test rig for heat-transfer studies of symmetrical impulse turbine 

blades . 
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1/32 -inch -outside -diameter 
Alundum tubing 


0.035- inch-diameter by 0.032- inch 
Inconel plug 


0.055- inch cutside- 
diam. Inconel tubing 




36 gage chrome l-alumel virei 


I 


C- 20285 
12.22.47 


Figure 4. - Radiation- type thermocouple probe 
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Figure 5. - Comparison of experimental blade temperature 
distribution and theoretical curve obtained by least 
squares. 0 = 6.78 cosh 4l t 68 (0^983-x) 
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1.0 



.4 .6 .8 

Exit Mach number, M Q 


Figure 6. - Recovery factor and aZ/PT for impulse 

turbine blades as functions of exit Mach number, 
Exit Reynolds number based on perimeter divided 
by 7 i. Gas properties taken at film temperature. 
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Reynolds number. Re 


Figure 7. - Correlation of hea t-t rans f e r data for air flowing 
past heated turbine blades. Pr, fig, kg based on film tem- 
perature. Characteristic blade dimension equal to perimeter 
divided by n. 
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Figure 8. - Heat-t rans t er data for three constant exit Reynolds 
numbers and variable exit Mach numbers. Characteristic blade 
dimension equal to perimeter divided by n. Gas properties 
taken at film temperature. 
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Figure 9. - Comparison of heat-transfer data for turbine-blade 
cascades, single cylinders, and streamline bodies. Reynolds 
and Nusselt numbers based on perimeter divided by n. Physical 
properties taken at film temperature. 
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F f g u re 10. - Comparison of theoretical and 
experimental blade temperature distribution 
at effective gas temperature of 1167° F„ 

• -9.73 cosh 43 ‘ 4 ( 1 • 25 ~* * 
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Figure II. - Coordinate system for radiation- 
probe temperature-distribution analysis. 
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